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Abstract: The requirement of t — b — t Yukawa coupling unification is common in sim- 
ple grand unified models based on the gauge group 50(10), and it also places a severe 
constraint on the expected spectrum of superpartners. For Yukawa-unified models with 
(j, > 0, the spectrum is characterized by three mass scales: i). first and second generation 
scalars in the multi-TeV range, ii). third generation scalars, \i and tua in the few-TeV range 
and Hi), gluinos in the ~ 350 — 500 GeV range with chargino masses around 100 — 160 
GeV. In such a scenario, gluino pair production should occur at large rates at the CERN 
LHC, followed by gluino three-body decays into neutralinos or charginos. Discovery of 
Yukawa-unified SUSY at the LHC should hence be possible with only 1 fb _1 of integrated 
luminosity, by tagging multi-jet events with 2-3 isolated leptons, without relying on miss- 
ing Et- A characteristic dilepton mass edge should easily be apparent above Standard 
Model background. Combining dileptons with 6-jets, along with the gluino pair produc- 
tion cross section information, should allow for gluino and neutralino mass reconstruction. 
A secondary corroborative signal should be visible at higher integrated luminosity in the 
X1X2 ~~ * 31 channel, and should exhibit the same dilepton mass edge as in the gluino 
cascade decay signal. 
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1. Introduction 



Models for new physics based on 50(10) grand unification and super symmetry are espe- 
cially compelling|Q]. ,50(10) grand unified theories (GUTs) allow for- in addition to gauge 
coupling unification- matter unification in the sense that all fields of a single generation of 
the Standard Model (SM) are embedded into the 16-dimensional spinor representation of 
,50(10). Since there are only 15 fields in a SM generation, the 16th element is occupied by 
a right-handed neutrino (RHN) field. The RHN is necessary for implementing the see-saw 
mechanism [Q] of generating neutrino masses. In addition, since the gauge group 50(10) 
is anomaly free, the ad-hoc but fortuitous triangle anomaly cancellation in the SM and 
in SU(5) GUTs is explained. The addition of supersymmetry to the model allows for a 
stabilization of the weak scale-GUT scale hierarchy, and is supported by the fact that the 
measured weak scale gauge couplings actually unify at Q = Mqut — 2 x 10 16 GeV under 
minimal supersymmetric standard model (MSSM) renormalization group evolution (RGE). 
Along with gauge and matter unification, the simplest 50(10) SUSY GUT models may 
also have Yukawa coupling unification, especially for the third generation. 

On the minus side, 50(10) SUSY GUT models in 4 dimensions are often unwieldy, 
due in part to the large Higgs representations which are needed to break the GUT gauge 
symmetry. In recent years, progress has been made in constructing SUSY GUT theories 
in 5 or 6 spacetime dimensions @. In these cases, the GUT symmetry can be broken by 
extra-dimensional compactification on for instance an orbifold. Such theories maintain 
many of the essential predictions of SUSY GUTs, while dispensing with GUT-breaking 
via large Higgs representations. In any case, if we consider that the MSSM arises from 
superstring theory as a 4-d effective theory below Mqut, then 6-7 extra dimensions will 
have to be compactified, and the description of string— >GUT— >MSSM is likely to be more 
complicated than any of the existing models. 

To avoid dealing with the unknown physics above the GUT scale, in this paper we will 
assume that an 50(10) SUSY GUT model breaks (either via the Higgs mechanism or via 
compactification of extra dimensions) to the SM gauge group at energy scale Q = Mqjjt- 
Below Mqut, we assume the MSSM is the correct effective field theory which describes 
nature. We will further assume that the superpotential above Mqut is of the form 

/ 3 ZV'ieV'ie^io H (1-1) 

so that the three third generation Yukawa couplings ft, ft and f T are unified at Mqut- 1 We 
will further make the reasonable assumption that the GUT scale soft SUSY breaking terms 
are constrained by the 50(10) gauge symmetry such that matter scalars have a common 
mass miQ, Higgs scalars have a common mass mio and there is a common trilinear soft 
parameter Aq. The bilinear soft term B can be traded as usual for tan/3, the ratio of Higgs 
field vevs, while the magnitude of the superpotential Higgs mass [i is determined in terms 

1 It is simple in this context to include as well the effect of a third generation neutrino Yukawa coupling 
/„; this effect has been shown to be small, although it can affect Yukawa coupling unification by a few per 
cent, if the right hand neutrino mass scale is within a few orders of magnitude of Mgc/tPI- 



- 1 - 



of M| via the electroweak symmetry breaking minimization conditions. Here, electroweak 
symmetry is broken radiatively (REWSB) due to the large top quark mass. 

In order to accomodate REWSB, it is well-known that the GUT scale Higgs soft 
masses must be split such that rrvjj n < Trijj , in order to fulfill the EWSB minimization 
conditions. Two possibilities™ full D-tevm splitting applied to all scalar masses (DT), or 
splitting applied only to Higgs soft terms (HS)- have been considered || |3j]. The latter HS 
case was found to yield a higher degree of Yukawa coupling unification when mie ~ 10 
TeV, so that this is the case we will consider here. We parametrize the Higgs splitting 
as m 2 H = m| =p 2Mp. The Higgs splitting can originate via a large near-GUT-scale 
threshold correction arising from the neutrino Yukawa coupling: see Ref. || for discussion. 
Thus, the Yukawa unified SUSY model is determined by the parameter space 

mi6, m 10 , Mp, m x / 2 , A , tan/5, sign(fi) (1.2) 

along with the top quark mass, which we take to be rat = 171 GeV. 2 

In previous work, the above parameter space was scanned over (via random scans |?], |3j 
and also by more efficient Markov Chain Monte Carlo (MCMC) scans Q) to search for 
Yukawa unified solutions. The quantity 

p max(f t , f b , f T ) 
mm(f t , fb, U) 

was examined, where solutions with R ~ 1 gave apparent Yukawa coupling unification. 
For superpotential Higgs mass parameter fj, > (as favored by (g — 2) M measurements), 
Yukawa unified solutions with R ~ 1 were found but only for special choices of GUT scale 
boundary conditions 0, g §, g 0, §]: 

A 2mi6, mio ~ l.2m 16 , (1.4) 

with mx/2 *C ?7ii6 and tan/? ~ 50. In fact, models with this sort of boundary condi- 
tions were derived even earlier in the context of inverted scalar mass hierarchy models 
(IMH) which attempt to reconcile suppression of flavor-changing and CP-violating pro- 
cesses with naturalness via multi-TeV first / second generation and sub- TeV scale third gen- 
eration scalarsjllj. The Yukawa-unified spectral solutions were thus found in Refs. |, g 



to occur with the above peculiar choice of boundary conditions along with mig ~ 3 — 20 
TeV and ~ 50 — 100 GeV. For accurate spectrum evaluation, we adopted the Isajet 
7.75 algorithm[12, 13], which includes 

• full two-loop RG evolution of both couplings and soft SUSY breaking (SSB) terms WM . 

• minimization of the RG-improved one-loop scalar potential at an optimized scale 
Msusy (we take Msjjsy = J^W L m l R ) > which accounts for leading two- loop terms [p^j, 

• implementation of complete 1-loop t, b and r threshold corrections |l6|, 17] and 

• implementation of complete 1-loop corrections for all sparticle masses [17]. 



2 Recent measurements from CDF and Do|| now find a central value for m t of 172.6 ± 1.4 GeV. 
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a hybrid approach of changing gauge and Yukawa coupling beta functions as different 
soft terms are integrated out of the effective theory pi 



A plot of Yukawa coupling evolution for case A of Table 1 of Ref. Q with m\% ~ 9 TeV 
is shown in Fig. [j]. Here, we clearly see the importance of the MSSM theshold corrections, 
which give rise to the discontinuities located around Q = Msusy ~ 3 TeV. 

Based on the above work[||, ||], the sparticle mass spectra from Yukawa-unified SUSY 
models are characterized by the following conditions: 

• first and second generation scalars have masses in the 3 — 20 TeV regime, 

• third generation scalars, /i and tua have masses in the 1 — 4 TeV regime (owing to 
the inverted scalar mass hierarchy), 

• gluinos g have mass in the 350 — 500 GeV range, 

• light charginos have mass ~ 100 — 160 GeV and 

• the lightest neutralino Xi is nearly pure bino with mass m^,o ~ 50 — 80 GeV. 

The presence of a bino-like Xi along with multi-TeV scalars gives rise to a neutralino 
dark matter (DM) relic abundance that is typically in the range ^^h 2 ~ 10 2 — 10 5 , i.e. 
far abovegH the WMAP measured|]| value VtcDNih 2 = O.lll+g^g (2a) by several orders 
of magnitude. 

A very compelling way out for solutions with m\§ ~ 5 — 20 TeV, is to not assume 
that the lightest neutralino Xi is the lightest SUSY particle (LSP). Instead, we assume 
that the axino a (the spin-^ element of the supermultiplet containing the axion) is in fact 



**-r 0.7 




le+00 le+02 le+04 le+06 le+08 le+10 le+12 le+14 le+16 



Q (GeV) 

Figure 1: Plot of ft, fb and f T evolution from the weak scale to the GUT scale for Point A of 
Table |. The large jumps around 3 TeV correspond to the MSSM threshold corrections. 
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the LSP[20|. The axino mass is very model-dependent, but can be anywhere in the keV- 
GeV range [^l], |2(J. In this case, the neutralino mass over-abundance is greatly reduced 
since Xi — * decay can occur with a lifetime of order 0.03 sec. This decay time is 
sufficiently short that late-time neutralino decay to axino in the early universe should not 
upset successful predictions of Big Bang Nucleosynthesis (BBN). The relic abundance then 
gets reduced by the ratio ma / m^o which can be of order 10~ 2 - 10" 5 . The a coming from 
Xi decay would actually constitute warm dark matter^]. However, axinos can also be 
produced thermally in the early universe as cold dark matter^, ^3|. Thus, in the axino 
LSP scenario, we would actually have dominantly thermally produced cold axino DM with 
a small admixture of warm axino DM arising from Xi ~~ > ^7 decays. In addition, there is 
also the possibility of a significant presence of axion CDM. 

In fact, in this class of Yukawa-unified solutions with miQ in the multi-TeV range, 
we also expect the gravitino G to lie in the multi-TeV range. The cosmological gravitino 
problem- wherein gravitinos produced thermally in the early universe suffer a late-time 
decay, thus destroying the successful predictions of BBN- can be avoided. For ~ 5 
TeV, the re-heat temperature Tr must be Tr ~ 10 5 GeV[24], thus creating tension with 
most viable mechanisms for baryogenesis [25]. However, for rriQ ~ 5 TeV, the re- heat 
bound is much higher: Tr ~ 10 8 — 10 9 GeV. This range of Tr is exactly what is needed 
for baryogenesis via non-thermal leptogenesis[26], wherein the heavy right-hand neutrino 
states are not produced thermally, but rather via inflaton decay. It also turns out to be 
the exact range needed to generate a dominantly cold axino DM universe. Thus, the whole 
scenario fits together to offer a consistent cosmological picture of BBN, baryogenesis and 
cold DM composed of axinos [27]! 

An alternative but equally compelling solution to the CDM problem in Yukawa-unified 
SUSY models occurs in cases where mie is as low as the ~ 3 TeV range. In this case, the 
Xi can remain as the LSP. Since tuiq is so low, Yukawa coupling unification only occurs at 
the R ~ 1.09 level. However, XiXi can now annihilate through the light Higgs h resonance 
at a sufficient rate to obtain the desired relic density with Xi remaining as the LSP||]. 

In this paper, we explore the consequences of Yukawa-unified SUSY models for sparticle 
detection at the LHC. We focus most of our attention on two cases presented in Table 1 of 
Ref. |||: 1. point A with m\Q ~ 9 GeV and an axino LSP, and 2. point D with m\Q ~ 3 TeV 
and a neutralino LSP. While our studies here focus on just two cases, we maintain that the 
qualitative features of all Yukawa-unified SUSY LHC signatures should be rather similar 
to these two cases. In fact, the collider phenomenology of these cases is rather similar 
between the two, since in the first case the neutralino decays to an axino far beyond the 
detector boundaries. Thus, in both cases the lightest neutralino Xi leads to missing Et at 
collider experiments. 

For the benefit of the reader, we present in Table [l] the two case studies we examine. We 
present the parameter space values, sparticle mass spectrum, and in addition the total tree- 
level LHC sparticle production cross section. We also list as percentages some contributing 
2^2 subprocess reactions. 

The remainder of this paper is organized as follows. In Sec. 0, we present some details 
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parameter 


Pt. A 


Pt. D 


mi 6 


9202.9 


2976.5 


mi/2 


62.5 


107.0 


A 


-19964.5 


-6060.3 


mio 


10966.1 


3787.9 


tan (3 


49.1 


49.05 


M D 


3504.4 


1020.8 


ft 


0.51 


0.48 


fb 


0.51 


0.47 


fr 


0.52 


0.52 


n 


4179.8 


331.0 


m~ g 


395.6 


387.7 


mu L 


HI OCT A 

9185.4 


2970.8 


m h 


zolo.l 


AO A £ 

434.5 


m b x 


0700 1 


o4y.o 


1TLe L 


Q1 ^1 Q 


zyoo.o 


m xt 


Izo.o 


1 7 
1UO. / 


128.6 


105.1 


X\ 

m A 


55.6 


52.6 


3273.6 


776.8 


m h 


125.4 


111.1 


a[fb] 


75579.1 


89666.1 


% m 


86.8 


80.5 


% (xfxi) 


8.8 


12.8 


% (Ml) 





1.1 



Table 1: Masses and parameters in GeV units for two cases studies points A and D of Ref. || 
using Isajet 7.75 with m t = 171.0 GeV. We also list the total tree level sparticle production cross 
section in fb at the LHC, plus the percent for several two-body final states. 



involving total sparticle production rates at the LHC and also sparticle branching fractions 
expected in Yukawa-unified SUSY models. In Sec. ||, we generate simulated LHC collider 
events associated with these two scenarios along with SM background (BG) rates. Since 
the collider signals are quite similar for both of points A and D, we mainly present figures 
for just point A. Our main findings include: the Yukawa-unified SUSY scenarios should be 
easily discoverable at LHC with about 1 fb -1 of integrated luminosity, without using missing 
Et cuts, by requiring events with a high isolated lepton multiplicity: ri£ > 2 or 3. In Sec. 
we present several distributions which can be used to gain information on sparticle masses, 
especially m g , m^o and m^o. The gluino mass m g can be extracted from total cross section 
results. In Sec. ||, we discuss a corroborating signal in the pp — > X1X2 ~* 3£+E™ lss channel. 
This channel comes from the soft component of signal events involving direct production 
of the lighter charginos and neutralinos. Our conclusions are presented in Sec. ||. 
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2. Cross sections and branching fractions for sparticles in Yukawa-unified 
models 



Given the characteristic spectrum of superpartners obtained in Yukawa-unified SUSY mod- 
els, it is useful to examine what sort of new physics signals we would expect at the LHC. 
Obviously, first /second generation squarks and sleptons in the multi-TeV mass range will 
essentially decouple from LHC physics. Gluinos- in the 350-500 GeV range- will be pro- 
duced in abundance via qq and gg fusion subprocesses. Charginos and neutralinos, being 
in the 100-160 GeV range, may also be produced with observable cross sections. 

As noted above, we list the tree-level total sparticle production cross sections obtained 
from Isajet for cases A and D in Table |l]. In case A, we find a (tot) ~ 8 x 10 4 fb, so that 
8000 sparticle pair events are expected at LHC with just 0.1 fb -1 of integrated luminosity. 
Of this total, 86.7% comes from gluino pair production, while 8.8% comes from X1X2 
production and 4.5% comes from X1X1 production. In case D, a(tot) ~ 9.6 x 10 4 fb, 
with 80.4% from gg production, 12.8% from X1X2 production, 6.4% from xfx.i production 
while top-squark pair production yields just 1.1% of the total. Given these production 
cross sections, we expect Yukawa-unified SUSY to yield primarily gg events at the LHC. 
Gluino pair production typically leads to events with hard jets, hard E™ 1SS and isolated 
leptons from the gluino cascade decays pS(|. We also expect a soft component coming from 
X1X1 an d xr~X2 production. While both these reactions lead to events with rather soft 
jets, leptons and E™ ss , the latter reaction can also yield clean trilepton events [29], which 
should be visible at LHC above SM backgrounds. 

For the case of gluino masses other than those listed in Table |l|, we show in Fig. || 
the total gluino pair production rate versus m g at the LHC at tree level (solid) and next- 
to-leading-order (NLO) using the Prospino program|30|]. The scale choice is taken to be 
Q = m g . We take rriq to be 3 TeV (blue) and 9 TeV (red). As can be seen, the results 
hardly vary between this range of squark masses. The tree level results agree well with 
Isajet, but the NLO results typically show an enhancement by a factor ~ 1.6. Thus, we 
expect Yukawa- unified SUSY models to yield pp — > gg X events at a 30-150 pb level at 
LHC. 

In Fig. g we show the total -ino pair production cross sections versus chargino mass 
m -± . While X1X2 an d X1X1 production dominate, and have rates around 10 3 — 10 4 fb over 
the range of interest, there exists a sub-dominant rate for xfxi anci a l so X1X2 production. 

Now that we see that Yukawa-unified SUSY will yield dominantly gluino pair pro- 
duction events at the LHC, we next turn to the gluino branching fractions in order to 
understand their event signatures. All sparticle branching fractions are calculated with 
Isajet 7.75. In Fig. we show various gluino branching fractions for points A and D. We 
see immediately that in both cases, BF(g — > bbx®) dominates at around 56%. This is fol- 
lowed by BF(g -► bb~Xi) at ~ 16%, and BF(g -> bixf) and BF(g -► tbxi) each at ~ 10 %- 
Decays to first and second generation quarks are much suppressed due to the large first 
and second generation squark masses. From these results, we expect gluino pair production 
events to be rich in 6-jets, E™ 1SS and occassional isolated leptons from the leptonic decays 
Xi — ► i^eXi an d X2 ~^ ^ + ^~x\, where I = e or (j,. 
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Figure 2: Plot of cr(pp — > ggX) in pb at y/s — 14 TeV versus m~ g . We use Prospino with scale 
choice Q = nig , and show LO (solid) and NLO (dashes) predictions in the vicinity of point A (red) 
and point D (blue) from Table |l|. 
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Figure 4: Plot of various sparticlc branching fractions taken from Isajet for points A and D from 
Table @. 

3. Gluino pair production signals at the LHC 



To examine collider signals from Yukawa-unified SUSY at the LHC in more detail, we 
generate 10 6 sparticle pair production events for points A and D, corresponding to 13 and 
11 fb _1 of integrated luminosities. We use Isajet 7.75|12] for the simulation of signal and 
background events at the LHC. A toy detector simulation is employed with calorimeter cell 
size Ar] x A(j> = 0.05 x 0.05 and —5 < r/ < 5. The HCAL energy resolution is taken to be 
80%/VE + 3% for \rj\ < 2.6 and FCAL is 100%/VE + 5% for \rj\ > 2.6. The ECAL energy 
resolution is assumed to be 3%/VE + 0.5%. We use a UAl-like jet finding algorithm with 
jet cone size R = 0.4 and require that Ex(jet) > 50 GeV and \rj(jet)\ < 3.0. Leptons are 
considered isolated if they have pr(e or fi) > 20 GeV and \rj\ < 2.5 with visible activity 
within a cone of AR < 0.2 of T,Ej? lls < 5 GeV. The strict isolation criterion helps reduce 
multi-lepton backgrounds from heavy quark (cc and bb) production. We also invoke a 
lepton identification efficiency of 75% for leptons with 20 GeV< pr(^) < 50 GeV, and 85% 
for leptons with pr(i) > 50 GeV. 

We identify a hadronic cluster with Et > 50 GeV and \rj(j)\ < 1.5 as a 6-jet if it 
contains a B hadron with pr(B) > 15 GeV and < 3 within a cone of AR < 0.5 

about the jet axis. We adopt a 6-jet tagging efficiency of 60%, and assume that light quark 
and gluon jets can be mis-tagged as 6-jets with a probability 1/150 for Et < 100 GeV, 



1/50 for E T > 250 GeV, with a linear interpolation for 100 GeV< E T < 250 GeV|3J. 



-8- 



process 


events 


cr (fb) 


Gl 


i i remiss 

G 1 + E T 


QCD (pt : 0.05 — 0.1 TeV) 


10 6 


2.6 x 10 


4.1 x 10 5 




QCD (p T : 0.1 - 0.2 TeV) 


10 6 


1.5 x 10 9 


1.4 x 10 7 




QCD (p T ■ 0.2 - 0.4 TeV) 


10 6 


7.3 x 10 7 


6.5 x 10 6 


2199 


k^LjD (pt : U.4 — l.U leVJ 


ID 


z.7 x 1U 


Z.5 X ID 


i i 

1157 


yUD (j?^ : 1 ~~ leVJ 


11) 


1 c . v 1 r\4 

l.o X ID 


lUoZ 


Zo 


TT7" /J i ■ « 

14/ — > tf£ + jets 


5 x 10 


3.9 x 10 


3850 


1275 


Z —> tt + jets 


5 x 10 


1.4 x 10 


1358 


652 


tt 


O X 10 


4.9 X 1U 


O.Z X 11) 


28 M 


WW, zz, wz 


5 x 10 5 


8.0 x 10 


197 


7 


lotal BG 


9.5 x 10 


2.76 x 10 


2.13 x 10 


8188 


Fomt A: 


10 


/.b x 10 


3.6 x 10 


8914 


S/B -> 






0.002 


1.09 


+ S (lfb- 1 ) -► 








68 


Point D: 


10 6 


9.0 x 10 4 


3.7 x 10 4 


10843 


S/B -» 






0.002 


1.32 


S/VS + B (lfb- 1 ) -» 








78 



Table 2: Events generated and cross sections (in fb) for various SM background and signal processes 
before and after cuts. The Gl' and E™ lss cuts are specified in the text. The W + jets and Z + jets 
background has been computed within the restriction pt{W, Z) > 100 GeV. 



In addition to signal, we have generated background events using Isajet for QCD jet 
production (jet-types include g, u, d, s, c and b quarks) over five pt ranges as shown 
in Table |2[ Additional jets are generated via parton showering from the initial and final 
state hard scattering subprocesses. We have also generated backgrounds in the W + jets, 
Z + jets, ti (with nit = 171 GeV) and WW, WZ, ZZ channels at the rates shown in 
Table ||. The W + jets and Z + jets backgrounds use exact matrix elements for one parton 
emission, but rely on the parton shower for subsequent emissions. 

First we require modest cuts: n(jets) > 4. Also, SUSY events are expected to spray 
large Et throughout the calorimeter, while QCD dijet events are expected to be typically 
back-to-back. Thus, we expect QCD background to be peaked at transverse sphericity 
St ~ 0, while SUSY events have larger values of St- 3 The actual St distribution 
for point A is shown in Fig. || (the St distribution for point D is practically the same). 
Motivated by this, we require St > 0.2 to reject QCD-like events. 

We plot the jet Et distributions of the four highest Et jets from Pt. A (in color) and 
the total SM background (gray histogram) in Fig. ^, ordered from highest to lowest Et, 

3 Here, St is the usual sphericity variable, restricted to the transverse plane, as is appropriate for hadron 
colliders. Sphericity matrix is given as 

\Y,PzPy 22 Py J 

from which St is defined as 2Ai/(Ai + A2), where An.,2 are the larger and smaller eigenvalues of S. 



-9- 



with jets labelled as jl — j4. The histograms are normalized to unity in order to clearly 
see the differences in distribution shapes. Again, the distributions for point D look very 
similar. We find that the highest Et jet distribution peaks around Et ~ 150 GeV with 
a long tail extending to higher Et values, while for the background it peaks at a lower 
value of Et ~ 100 GeV. Jet 2 and jet 3 have peak distributions around 100 GeV both for 
the signal and backgrounds, while the jet 4 distribution backs up against the minimum jet 
Et requirement that ET(jet) > 50 GeV. Thus, at little cost to signal but with large BG 
rejection, we require Et{j1) > 100 GeV. 



The collection of cuts so far is dubbed Cl'||32j: 



CI' cuts: 



n(jets) > 4, 
E T (jl,j2J3,j4) > 100, 50,50,50 GeV, 
S T > 0.2. 



(3.2) 
(3.3) 
(3.4) 



The classic signature for SUSY collider events is the presence of jets plus large E™ ss p3| . 
In Fig. [?], we show the expected distribution of .E™ 1SS from points A and D, along with SM 
BG. We do see that signal becomes comparable to BG around £?p lss ~ 150 GeV. We list 
cross sections from the two signal cases plus SM backgrounds in Table ||| after cuts Cl' plus 
£vmiss y 150 GeV. While signal S is somewhat higher than the summed BG B, the signal 
and BG rates are rather comparable in this case: S/B = 1.09 for pt. A while S/B = 1.32 
for pt. D. 

Even so, it has been noted in Ref. [^] that E™ ss may be a difficult variable to reliably 
construct during the early stages of LHC running. The reason is that missing transverse 
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Figure 5: Plot of distribution in transverse sphericity St for events with cuts Cl' from benchmark 
point A and the summed SM background; point D leads to practically the same distribution. 
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energy can arise not only from the presence of weakly interacting neutral particles such as 
neutrinos or the lightest neutralinos, but also from a variety of other sources, including: 

• energy loss from cracks and un-instrumented regions of the detector, 

• energy loss from dead cells, 

• hot cells in the calorimeter that report an energy deposition even if there is not one, 

• mis- measurement in the electromagnetic calorimeters, hadronic calorimeters or muon 
detectors, 

• real missing transverse energy produced in jets due to semi-leptonic decays of heavy 
flavors, 

• muons and 

• the presence of mis- identified cosmic rays in events. 

Thus, in order to have a solid grasp of expected E™ 1SS from SM background processes, it 
will be necessary to have detailed knowledge of the complete detector performance. As 
experience from the Tevatron suggests, this complicated task may well take some time to 




200 300 



400 500 

E T (/.) (GeV) 



500 



Figure 6: Plot of jet Et distributions for events with > 4 jets after requiring just St > 0.2, from 
benchmark point A; distributions for point D are nearly the same. 
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complete. The same may also be true at the LHC, as many SM processes will have to 
be scrutinized first in order to properly calibrate the detector [34]. For this reason, SUSY 
searches using E™ ss as a crucial requirement may well take rather longer than a year to 
provide reliable results. 

For this reason, Ref. [32| advocated to look for SUSY signal events by searching for 
a high multiplicity of detected objects, rather than inferred undetected objects, such as 
E™ 1SS . In this vein, we show in Fig. [8| the jet multiplicity from SUSY signal (Pts. A and D) 
along with SM BG after cuts CI', i.e with no E™ 1SS cut. We see that at low jet multiplicity, 
SM BG dominates the SUSY signal. However, signal/background increases with n(jets) 
until at n(jets) ~ 15 finally signal overtakes BG in raw rate. 

One can do better in detected 6-jet multiplicity, n&, as shown in Fig. ||. Since each 
gluino is expected to decay to two 6-jets, we expect a high multiplicity in signal. In this 
case, BG dominates signal at low n&, but signal overtakes BG around ~ 4. 

The isolated lepton multiplicity ni is shown in Fig. 10 for signal and SM BG after cuts 
CI'. In this case, isolated leptons should be relatively common in gluino cascade decays. 
We see that signal exceeds BG already at ri£ = 2, and far exceeds BG at ne = 3: at ri£ = 3, 
we have 35 (42) fb of signal cross section after cuts for point A (D) and no BG. In fact, 
high isolated lepton multiplicity was advocated in Ref. [j32j l in lieu of an E™ 1SS cut to search 
for SUSY with integrated luminosities of around 1 fb -1 at LHC. 

We also point out here that gg production can lead to large rates for same-sign (SS) 
isolated dilepton production] 35 1, while SM BG for this topology is expected to be small. 



We plot in Fig. 11 the rate of events from signal and SM BG for cuts CI' plus a pair 
of isolated SS dileptons, versus jet multiplicity. While BG is large at low n(jets), signal 
emerges from and dominates BG at higher jet multiplicities. 




Figure 7: Plot of missing Et for events with > 4 jets after cuts CI', from benchmark points A 
(full red line) and D (dashed blue line). 
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Figure 8: Plot of jet multiplicity from benchmark points A (full red line) and D (dashed blue line) 
after cuts CI' along with SM backgrounds. 
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Figure 9: Plot of 6-jet multiplicity from benchmark points A (full red line) and D (dashed blue 
line) after cuts CI' along with SM backgrounds. 
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Figure 10: Plot of isolated lepton multiplicity from benchmark points A (full red line) and D 
(dashed blue line) after cuts CV along with SM backgrounds. 



Cuts S > 0.2 + 2 SS dileptons 




Figure 11: Plot of jet multiplicity in events with isolated SS dileptons from benchmark points A 
(full red line) and D (dashed blue line) after cut St > 0.2 along with SM backgrounds. 
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4. Sparticle masses from gluino pair production 



There exists good prospects for sparticle mass measurements in Yukawa-unified SUSY 
models at the LHC. One reason is that sparticle pair production is dominated by a single 
reaction: gluino pair production. The other propitious circumstance is that the mass 
difference m^o — m^o is highly favored to be bounded by Mz- This means that x\ decays 
dominantly into three body modes such as X2 ~^ Xi^ with a significant branching fraction, 
while the so-called "spoiler decay modes" X2 ~^ Xi% an d X.2 ~^ Xih are kinematically 
closed. The three body decay mode is important in that it yields a continuous distribution 
in m(££) which is bounded by m^o — m^o: this kinematic edge can serve as the starting 
point for sparticle mass reconstruction in cascade decay events p^, |29|| . 

As an example, we require cuts set Cl' plus the presence of a pair of same-flavor 
opposite-sign (SF/OS) isolated leptons. The resulting dilepton invariant mass distributions 
for points A and D are shown in Fig. |l2[ Furthermore, in Fig. |l^, we plot the different-flavor 
subtracted distributions: da/dm(£ + £~) — da/dm(£ + £~'), which allow for a subtraction of 
e + fi~ and e~ fj, + pairs from processes like chargino pair production in cascade decay events. 
A clear peak at m(££) = Mz is seen in the BG distribution. This comes mainly from 
QCD jet production events, since Isajet includes W and Z radiation in its parton shower 
algorithm (in the effective W approximation) . The signal displays a histogram easily visible 
above SM BG with a distinct cut-off at m^o — m^o = 73 (52.5) GeV for point A (D). Isajet 
contains the exact decay matrix elements in 3-body decay processes, and in these cases we 
see distributions that differ from pure phase space, and yield distributions skewed to higher 
m(££) values. This actually shows the influence of the virtual Z in the decay diagrams, 
since the decay distribution is dominated by Z* exchange. The closer m^o — m^o gets to 
Mz, the more the Z-boson propagator pulls the dilepton mass distribution towards Mz|j37||. 
The dilepton mass edge should be measureable to a precision of ~ 50 MeV according to 
Ref. ||]. 

For Yukawa-unified SUSY models, the branching fraction BF(g — > 66x2) dominates 
at around 56%. If one can identify events with a clean \\ ~^ ^Xi decay, then one might 
also try to extract the invariant mass of the associated two 6-jets coming from the gluino 
decay, which should have a kinematic upper edge at m g — m^o ~ 267 (283) GeV for point 
A (D). A second, less pronounced endpoint is expected at m g — m^o ~ 340 (335) GeV due 
to g —> bbxi decays which have ~ 16% branching ratio. A third endpoint can also occur 
from X2 ~ * Xibb decay where m^o — m^o = 73 (52.5) GeV, respectively. 

The high multiplicity of 6-jets (typically two from each gluino decay) however poses a 
serious combinatorics problem in extracting the bb invariant-mass distribution. In a first 
attempt, we required at least two tagged 6-jets along with cuts Cl' and SF/OS dileptons, 
and plotted the minimum invariant mass of all the 6-jets in the event. The resulting 
distributions for points A and D peaked around m(bb) ~ 100 GeV, with a distribution 
tail extending well beyond the above-mentioned kinematic endpoints. It was clear that we 
were frequently picking up wrong 6-jet pairs, with either each b originating from a different 
gluino, or one or both b's originating from a x§ -> Xi^b decay. A more sophisticated 
procedure to pair the correct 6-jets is required for the bb invariant-mass distribution. 
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Figure 12: SF/OS dilepton invariant mass distribution after cuts CI' from benchmark points A 
(full red line) and D (dashed blue line) along with SM backgrounds. 
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Figure 14: Plot of m(Xi^) from benchmark points A and D along with SM backgrounds in events 
with cuts CI' plus > 4 &-jets and minimizing Am(Xi — X2); see text for details. 



Parton-level Monte Carlo simulations revealed that the two hardest (highest Er) b- 
jets almost always originated from different gluinos. Thus, we require events with at least 
four tagged 6-jets (along with cuts CI') and combine the hardest 6-jet with either the 3rd 
or 4th hardest 6-jet, creating an object X\(bb). Moreover, we combine the 2nd hardest 
6-jet with the 4th or 3rd hardest 6-jet, creating an object X 2 (bb). Next, we calculate 
AmpTi - X 2 ) = \m(Xi) - m(X 2 )\/(rn(Xi) + m(X 2 )). We select the set of bb clusters 
which has the minimum value of Am(Xi — X 2 ), and plot the invariant mass of both 
the clusters. This procedure produces a sharp kinematic edge in m(bb) in parton level 
simulations. 4 The resulting distribution from Isajet is shown in Fig. [l4| for points A and D. 
The distribution peaks at a higher value of m(Xi) (i = 1, 2), and is largely bounded by the 
kinematic endpoints, although a tail still extends to high m{Xj). Part of the high m(Xf) 
tail is due to the presence of g — > yf{bb decays, which have a higher kinematic endpoint 
than the g — ► x 2 bb decays, and of x 2 ~^ Xibb decays. In addition, there is a non-negligible 
background contribution, indicated by the gray histogram. 

We can do much better, albeit with reduced statistics, by requiring in addition the 
presence of a pair of SF/OS dileptons. Applying the same procedure as described above, 



we arrive at the distribution shown in Fig. 15. In this case, the SM BG is greatly reduced, 
and two mass edges begin to appear. 

It should also be possible to combine the invariant mass of the SF/OS dilepton pair 



4 We have also tried other methods such as picking the X\ — X2 pair with maximum A0(Xi — X2), 
maximum /S.R(Xi — X%), minimum 8pr{X\ — X2) and minimum average invariant mass avg(m(Xi), m(Xa)). 
We also tried to separate the resulting 6 jets into hemispheres. In the end, the best amount of correct 
assignment was achieved with the choice of the Am(Xi — X2)min- 
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Figure 15: Same as Fig. [l4| but requiring in addition a pair of SF/OS leptons. 



with a bb pair. Requiring cuts CI' plus > 2 6-jets plus a pair of SF/OS dileptons (with 
m(££) < m^o — m^o), we reconstruct m{bb££). The result is shown in Fig. [l(]. While the 
distribution peaks at m(££bb) ~ 300 GeV, a kinematic edge at rrig — m^o ~ 340 GeV is also 
visible (along with a mis-identification tail extending to higher invarant masses). 

We try to do better, again with a loss of statistics, by requiring > 4 6-jets instead of 
> 2, and combining the bb clusters into objects X\ and X2 as described above. We again 
choose the set of clusters which give the minimum of Am(Xi — X2) and combine each 
of these clusters with the ii pair. We take the minimum of the two m(Xi££) values, and 
plot the distribution in Fig. 17. In this case, the SM BG is even more reduced, and the 
771,5 — m v° mass edge seems somewhat more apparent. 

According to [38], measurements of hadronic mass edges can be made with a precision 
of roughly 10%. Nevertheless, from the kinematic distributions discussed above we can 
only determine mass differences. There is still not enough information to extract absolute 
masses, i.e. each of m,g, m^o and 777,^0. However, it is pointed out in Ref. |39| that in cases 
(such as the focus point region of minimal supergravity) where sparticle pair production 
occurs nearly purely from gg production, and when the dominant g branching fractions 
are known (from a combination of theory and experiment), then the total gg production 
cross section after cuts allows for an absolute measurement of 777^ to about an 8% accuracy. 
These conditions should apply to our Yukawa-unified SUSY cases, if we assume the ~ 56% 
branching fraction for g — > 66x2 decay (from theory). The study of Ref. [39| required that 
one fulfill the cuts C2 which gave robust gluino pair production signal along with small 
SM backgrounds: 
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Figure 16: Plot of m(bb£ + £ ) m in from points A and D along with SM backgrounds. 
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Figure 17: Plot of m(Xifll + £ ) m in from points A and D, minimizing Am(Xi — X2) as explained 
in the text, along with SM backgrounds. 
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C2 cuts: 



£^ iss > (max(100 GeV, 0.2M e// ), (4.1) 

n(jets) > 7, (4.2) 

n(b-jets)>2, (4.3) 

E T (jl,j2 - j7) > 100, 50 GeV, (4.4) 

A T > 1400 GeV, (4.5) 

S T > 0.2 , (4.6) 

where At is the augmented effective mass At = E^ 1 ^ + Ylieptons + Yljets ^T- In this 



case, the summed SM background was about 1.6 fb, while signal rate for Point A (D) 
is 57.3 (66.2) fb. The total cross section after cuts varies strongly with m g , allowing an 
extraction of m g to about 8% for 100 fb" 1 integrated luminosity, after factoring in QCD 
and branching fraction uncertainties in the total rate. Once an absolute value of m g is 
known, then m^o and m^o can be extracted to about 10% accuracy from the invariant 
mass edge information. 



5. Trilepton signal from X1X2 production 

While the signal from gluino pair production at the LHC from Yukawa-unified SUSY models 
will be very robust, it will be useful to have a confirming SUSY signal in an alternative 
channel. From Fig. ||| we see that there also exists substantial cross sections for xtxii 
xfXi and X1X.2 production. The xf -> Xiff an d xl ~> Xiff decays (here / stands 
for any of the SM fermions) are dominated by W and Z exchange, respectively, so that 
in this case the branching fractions BF(x^ — > Xiff') are similar to BF(W^ — > //') and 
BF(x% - xiff) is similar to Z - //". 

The xtxl XiQQ f + Xi process will be difficult to observe at LHC since the final state 
jets and £vp lss will be relatively soft, and likely buried under SM background. Likewise, 
the XiXi ~ * Xi^eXi signal will be buried under a huge BG from W — ► ivg production. 
The XiXi production reaction will also be difficult to see at LHC. The purely hadronic 
final state will likely be buried under QCD and Z + jets BG, while the lepton plus jets 
final state will be buried under W + jets BG. The dilepton final state will be difficult to 
extract from W + W~ and ti production. 

The remaining reaction, x^X^ production, yields a trilepton final state from X\ ~^ 
Xylvi an d x\ ~^ Xi^ decays which in many cases is observable above SM BG. The LHC 
reach for xfx^ ~^ 3£ + E™ 1SS production was mapped out in Ref. [4C], and the reach was 



extended into the hyperbolic branch/focus point (HB/FP) region in Ref. El. The method 
was to use the cut set SC2 from Ref. |37| but as applied to the LHC. For the clean trilepton 
signal from X1X2 ~~ * 3£ + -E™ 1SS production, we require: 

• three isolated leptons with pr{£) > 20 GeV and \r](\ < 2.5, 

• SF/OS dilepton mass 20 GeV < m(£ + £-) < 81 GeV, to avoid BG from photon and 
Z poles in the 2 — ► 4 process qq' — ► lit'vi , 
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Table 3: Clean trilepton signal after cuts listed in the text. 



• a transverse mass veto 60 GeV < Mt(£, 5™ 1SS ) < 85 GeV to reject on-shell W 
contributions, 

• Eft lss > 25 GeV and, 

• veto events with n(jets) > 1. 

The resulting BG levels and signal rates for points A and D are listed in Table |3|. 
The 2^2 processes are calculated with Isajet, while the 2 — ► 4 processes are calculated 
at parton level using Madgraphl[42|. The combination of hard lepton px cuts and the 



requirement that n(jets) = leaves us with no 2 — > 2 background, while the parton level 
2 — ► 4 BG remains at 0.7 fb. Here, we see that signal from the two Yukawa-unified points 
well exceeds background. 

In the clean 3£ channel, since two of the leptons ought to come from %2 — * X?^ decay, 
they should display a confirmatory dilepton mass edge at m^o — m^o as is evident in the 
gluino pair production events, where the dileptons are accompanied by high jet multiplicity. 
The distribution in m(£ + l~) is shown in Fig. |l^. Event rates are seen to be lower than 
those from gg production. Integrated luminosity needed for a discovery with 5a significance 
would be 2.83 fb" 1 for pt. A and 1.5 fb" 1 for pt. D 5 . 



6. Summary and conclusions 

Simple SUSY grand unified models based on the gauge group 5O(10) may have t — b — r 
Yukawa coupling unification in addition to gauge group and matter unification. By assum- 
ing the MSSM is the effective field theory valid below Mgut, we can, starting with weak 
scale fermion masses as boundary conditions, check whether or not these third generation 
Yukawa couplings actually unify. The calculation depends sensitively on the entire SUSY 
particle mass spectrum, mainly through radiative corrections to the 6, t and r masses. It 

5 Significance is defined as S/y/(S + B) 
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in the clean trilepton channel from points A and D along with SM 



was found in previous works that t — b — r Yukawa coupling unification can occur, but 
only for very restrictive soft SUSY breaking parameter boundary conditions valid at the 
GUT scale, leading to a radiatively induced inverted mass hierarchy amongst the sfermion 
masses. While squarks and sleptons are expected to be quite heavy, gluinos, winos and 
binos are expected to be quite light, and will be produced at large rates at the CERN LHC. 

We expect LHC collider events from Yukawa-unified SUSY models to be dominated 
by gluino pair production at rates of (30 — 150) x 10 3 fb. The gs decay via 3-body modes 
into bbx®, bb~Xi and tbxf, followed by leptonic or hadronic 3-body decays of the x\ and 
xf ■ A detailed simulation of signal and SM BG processes shows that signal should be 
easily visible above SM BG in the > 4 jets plus > 3£ channel, even without using the £vp lss 
variable, with about 1 fb -1 of integrated luminosity. 

If Yukawa-unified signals from gg production are present, then at higher integrated 
luminosities, mass edges in the m(£ + l~~), m{bb) and m(bb£ + £~) channels along with total 
cross section rates (which depend sensitively on the value of nig) should allow for sparticle 
mass reconstruction of m g , m^o and m^o to 0(10%) accuracy for ~ 100 fb -1 of integrated 
luminosity. The gluino pair production signal can be corroborated by another signal in 
the clean trilepton channel from X1X2 ~~ * 3£ + E™ ss , which should also be visible at higher 
integrated luminosities. Thus, based on the study presented here, we expect LHC to either 
discover or rule out t — b — r Yukawa-unified SUSY models within the first year or two of 
operation. 
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